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ABSTRACT 

One of the largest sources of uncertainty in the computation of globular cluster 
isochrones and hence in the age determination of globular clusters is the lack of 
a rigorous description of convection. Therefore, we calibrated the superadiabatic 
temperature gradient in the envelope of metal-poor low-mass stars according to the 
results from a new grid of 2D hydrodynamical models, which cover the Main Sequence 
and the lower Red Giant Branch of globular cluster stars. In practice, we still use for 
computing the evolutionary stellar models the traditional mixing length formalism, 
but we fix the mixing length parameter a in order to reproduce the run of the entropy 
of the deeper adiabatic region of the stellar envelopes with effective temperature and 
gravity as obtained from the hydro-models. The detailed behaviour of the calibrated 
a depends in a non-trivial way on the effective temperature, gravity and metallicity of 
the star. Nevertheless, the resulting isochrones for the relevant age range of galactic 
globular clusters have only small differences with respect to isochrones computed 
adopting a constant solar calibrated value of the mixing length. Accordingly, the age 
of globular clusters is reduced by 0.2 Gyr at most. 



Subject headings: convection — stars: evolution — stars: Population II — globular 
clusters: general 



1. Introduction 



One of the most important unsolved problems of stellar evolution is the determination of the 
temperature gradient in the superadiabatic regions at the top of the convective envelopes of cool 
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stars, which strongly affects the effective temperature (T c g) of these objects. The mixing length 
theory (MLT - Bohm-Vitense 1958) is widely used for deriving this gradient. It contains a number 
of free parameters, among them a, the ratio of the mixing length to the pressure scale height, 
which provides the scale length of the convection. There are different versions of the MLT, each 
one assuming different values for these parameters. As demonstrated by Pedersen et al. (1990), the 
T e fj values obtained from the different formalisms are equivalent, provided that a suitable value of 
a is selected (see also Gough & Weiss 1976). This means that the MLT results concerning stellar 
structure models depend only on one free parameter, namely a, and its absolute value depends on 
the selected MLT formalism. Once the formalism is fixed, a is usually calibrated by reproducing 
the solar T e g, and this solar-calibrated a is then used for computing models of stars very different 
from the Sun (e.g., metal-poor Red Giant Branch and Main Sequence stars of various masses). 
However, in principle there is no compelling reason that a should be the same for the Sun and 
different kinds of stars. 

More recently, Canuto & Mazzitelli (1991 - CM) proposed a new formalism for the treatment 
of the superadiabatic convection; they take into account the full turbulent energy spectrum and set 
the convective scale length equal to the geometrical depth from the top of the convective region. 
Comparisons between MLT (solar-calibrated a) and CM stellar models show that isochrones 
computed with the CM formalism cannot be reproduced by the MLT with any constant value of 
a (Mazzitelli et al. 1995). 

The problem of determining accurate effective temperatures for cool stars affects the globular 
cluster (GC) age determination and, in turn, the estimated age of the universe. Large variations of 
a alter the derived stellar T e g and colours and change the shape of GC isochrones; as demonstrated 
by Chaboyer (1995 - see also Chaboyer et al. 1998) the change of the isochrone shape can even 
modify the luminosity of the Turn Off (TO - bluest point along a given isochrone). Moreover, 
when comparing CM and MLT isochrones one finds that for the most metal-poor isochrones, in 
the relevant range of ages of the metal-poor galactic GCs, the TO luminosity obtained from the 
CM isochrones differs appreciably from the case of the MLT (Mazzitelli et al. 1995). Since the TO 
brightness is the main age indicator for GCs, uncertainties in the convection treatment can affect 
the estimated GC ages (by «1 Gyr or even more). 

In principle, one could try to constrain the convective efficiency by comparing theoretical 
isochrones with regions of the observed Colour-Magnitude diagrams of GCs whose colours are 
unaffected by the cluster age, but the still existing large uncertainties in the colour-transformations 
(see, e.g., Weiss & Salaris 1998) do not permit to safely follow this approach. An empirical 
constraint is given by the T e g of the upper Red Giant Branches of a sample of GCs as derived 
by Frogel et al. (1981). MLT models computed with a solar calibrated a (see, e.g., Vandenberg 
et al. 1996, Salaris & Weiss 1998) appear to be in agreement with these empirical data, also if 
a precise error bar on the a value calibrated in this way is hard to establish, being probably in 
the range ±0.1-0.3 (see Vandenberg et al. 1996 and references therein). However, this 'empirical' 
calibration in principle does not constrain the convection along the Main Sequence and the lower 
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Red Giant Branch of GCs. 

Both the MLT and CM formalisms assume a simplified, time-independent, one-dimensional, 
local treatment of a typically non-stationary, multi-dimensional, and non-local phenomenon 
which the stellar convection actually is. The final solution to the problem of the superadiabatic 
convection in stellar envelopes has to come from the computation of realistic multidimensional 
radiation-hydrodynamics (RHD) simulations covering the range of effective temperatures, gravities, 
and compositions typical of stars with convective envelopes. First attempts to include in stellar 
models the results from rather crude 2-dimensional (2D) and 3-dimensional (3D) hydrodynamical 
simulations date back to the works by Deupree & Varner (1980) and Lydon et al. (1992, 1993a,b). 

In this Letter we discuss the first application to the computation of GC isochrones of new 
results obtained from a grid of detailed 2D radiation-hydrodynamics models including realistic 
microphysics and a detailed treatment of the radiative transport. A further important feature of 
these models is that they span a wide range in metallicity ([M/H]=— 2.0 to [M/H]=0.0; here we 
adopt the usual spectroscopic notation [M/B]=log(M/ H) star -log(M/ H) Q , where M is the global 
metal abundance) and cover the Main Sequence (MS) and lower Red Giant Branch (RGB) region 
of GC colour- magnitude diagrams. 

2. Hydrodynamical models and calibration of the efficiency of convection 

The full grid of RHD models is described in detail elsewhere (Ludwig et al. 1998, Freytag et 
al. 1998a,b), and a comprehensive discussion about the numerical and physical assumptions of the 
RHD simulations can be found in Ludwig et al. (1994). Here we just recall the main features of 
the models. Each 2D model describes the atmosphere and upper layers of a star with a convective 
envelope. It is obtained by solving the time-dependent, non-linear equations of hydrodynamics 
for a stratified compressible fluid. The calculations take into account a detailed treatment of the 
equation of state and of the multi-dimensional, non-local, radiative transfer (for more details 
see Ludwig et al. 1994). Similar to classical model atmosphere calculations, the hydrodynamical 
models are fully determined by specifying T c g, acceleration of gravity g, and chemical composition, 
and they lie in the range 4300 K<T cff <7100 K, 2.54<log(#)<4.74, -2.0<[M/H]<0.0. 

From this grid of models one can extract the entropy of the deeper, adiabatic convective 
layers ( function of T e gf, g, and [M/H] (see Ludwig et al. 1998, Freytag et al. 1998a, b). 

Once this relation is implemented in a stellar evolution code, it completely fixes the T e g of the 
star as determined from the solution of the stellar structure equations. 

A way for implementing easily this dependence of s env on T e g- and g into an evolutionary 
code makes use of the MLT formalism. As explained in detail by Ludwig et al. (1998), for each 
fixed metallicity one can compute a grid of hydrostatic one-dimensional stellar envelope models 
based on the MLT, covering the same range of g and T e g spanned by the RHD computations, 
and using the same input physics. By employing as surface boundary condition the T(r) relation 
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derived from the hydro- models, one can calibrate an effective a (a c g) that is able to reproduce 
the 

Senv _ ^efT relation obtained from the RHD computations. In this way one can derive a function 
a efi = f(T e f[ ,g) a t each metallicity, that is easy to use for computing stellar evolutionary models. 
The estimated error on the values of a e g derived by means of this procedure is equal to ±0.05 
(Ludwig et al. 1998). 

Ludwig et al. (1998) discuss the comparison of their RHD model for the solar envelope with 
the results from helioseismology. They show that the entropy at the bottom of the superadiabatic 
region as derived from helioseismology would imply a value of a e ff for the Sun slightly higher, by 
«0.10±0.05, than the value deduced from their RHD models. This small discrepancy is explained 
by comparing their result with the outcome of similar 3D simulations and by the examination 
of the opacities used in their models. For the Sun the 3D models predict an increase of a by 
«0.07±0.02 with respect to the 2D ones, and the ATLAS6 opacities (Kurucz 1979) used in the 
models do not consider the contribution of the molecules. The effect of including this contribution 
to the opacity would further change a by wO.l. The combination of these two effects explains the 
small discrepancy between the adopted solar RHD model and the results from helioseismology. 
Another important result derived from the RHD models is that the effect of the envelope He 
abundance on the derived a e g values is basically negligible. 

The calibration of a c g for metal poor stars has been performed by Freytag et al. (1998a,b), 
and we have used their results for computing isochrones with typical GC metallicities [M/H]=— 2.0 
and [M/H]=— 1.0 (scaled solar metal distribution), Y=0.23, age t ranging from 9 to 14Gyr, using 
the code described in Salaris et al. (1997). We have employed the same T(r) relation and the 
same MLT formalism (Bohm-Vitense 1958) used in the calibration of a e g. 

We have computed a first set of isochrones using the ATLAS6 low temperature opacities. The 
only source of possible inconsistencies with the RHD models was in this case the equation of state 
(EOS) employed in the evolutionary calculations (see Salaris et al. 1997), which is not the same as 
in the RHD computations. Nevertheless, we have verified that it does not modify appreciably a e g 
as derived from the RHD calibration. For the sake of comparison, we have computed isochrones 
for the same age and metallicities, but using a constant, solar calibrated value of a (a©). Since 
the ATLAS6 data do not include the contribution of the molecules to the opacity of the stellar 
matter, we have repeated the same evolutionary computations previously described (with a c g and 
cxq), using this time the updated Alexander & Ferguson (1994) low temperature opacities, which 
include the molecular contribution. When using these opacity data, we found a small effect only 
on the zero point of the a e s=f(T c f{,g) relation. To fix the ideas, the solar calibration with the 
evolutionary code yields in this case a=1.69, while from the RHD models one gets a=1.59 for the 
Sun, that means a deviation by a factor 1.06. This small correction factor (which makes consistent 
the entropy of the adiabatic layers as derived from RHD models and from helioseismology) for the 
a cfr=f(7eff >fiO relation has therefore been taken into account in the evolutionary calculations. 

At this point we have verified that the differences among isochrones computed with a e s and 
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a are exactly the same in the case of models computed with ATLAS6 or Alexander & Ferguson 
(1994) opacities. Since the latter data are a more realistic evaluation of the opacity of stellar 
matter, in the following section we will discuss the isochrones computed using the Alexander & 
Ferguson (1994) opacities. 

Before concluding this section we would like to stress the fact that the derived calibration 
of a c g is only intended to reproduce the function s env (T e s,g) of the RHD models. The detailed 
temperature profile and convective velocities of the superadiabatic layers are not represented 
adequately by the MLT with a c g, but our main concern here is only the determination of reliable 
effective temperatures for cool stars. 



3. Results and discussion 

Representative isochrones (9 and 13 Gyr) for the two considered metallicities, computed with 
a c g and a Q are displayed in Figs. 1 and 2, upper panels. The most striking feature is the close 
resemblance between the two sets of isochrones. The MS loci are coincident, the T e g values of 
the TO points are very similar (notice the linear scale for the T e g axis), the biggest difference 
being equal to «45 K for the 9 Gyr most metal-poor isochrone (an age possibly too young for the 
metal-poor galactic GC population, see e.g. Salaris & Weiss 1998). Along the RGB the isochrones 
computed by using a c g are systematically hotter by only «50K for [M/H]=— 1.0 and «40K 
for [M/H]=— 2.0. To explain this behaviour it is useful to study the run of a c s with respect to 
log(L/L Q ) along the same isochrones, as shown in the lower panels of the same figures. 

The differences between a e g and a© along the lower MS are hardly relevant, since the T e g of 
these stars is insensitive to the choice of a (the entropy jump from the photosphere to the deep 
adiabatically stratified layers is small anyway), while around the TO they depend on the age of 
the isochrones. In general, in the youngest, most metal-poor isochrones a e g shows the largest 
difference with respect to a , but since stars in these phases are quite hot (T e g «7000 K) and 
their convection zones are relatively shallow, the sensitivity of T e s to a is not very large. Along 
the RGB, where the T e g- of stellar models is most sensitive to a because of deeper superadiabatic 
regions, a e ff is systematically higher than a Q by 0.10-0.15 for both metallicities. This difference 
causes a systematic shift by ^50 K toward higher T e ff with respect to the case of a & , a quantity 
marginally significant since the error by ±0.05 on a c g translates into an error by «±15-20K on 
the RGB T e ff). Qualitatively, the behaviour of isochrones computed using a e g looks similar, for 
certain features, to the results of the CM formalism; we are referring here to the fact that the TO 
is cooler (but only for the youngest more metal-poor isochrones) than for the models computed 
with otQ . But the differences we find are smaller than the predictions of the CM formalism. 
Moreover, the RGB location in the a e ff isochrones is only slightly hotter than in the «q ones, 
while in the case of CM models the RGB is cooler at low metallicity and progressively hotter for 
increasing metallicities. 
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At this point, let's turn our attention to the GC age indicators that one can extract from 
the isochrones. The TO brightness is the most solid one; once the distance is fixed (e.g., from 
the Horizontal Branch luminosity or by means of the subdwarf fitting technique), the comparison 
between theoretical and observed TO gives directly the cluster age. The TO colour is also a 
possible age indicator once the reddening is known, but the present uncertainties on the colour 
transformations do not favour this method for deriving absolute ages, also if the isochrone T e g and 
the GC reddenings are determined with high accuracy. A third possibility is to use the reddening- 
and distance modulus- independent quantity A(B — V), that is the difference in (B-V) between 
TO and base of the RGB, as defined by Vandenberg et al. (1990). Again, the uncertainties on the 
colour-transformations prevents from using the absolute value of A(B — V) for deriving absolute 
GC ages, but the differential use of this quantity is a solid and widely employed indicator of age 
differences (see, e.g., Vandenberg et al. 1990, Salaris & Weiss 1997), and it is weakly dependent 
on [M/H]. 

In Figure 3 (upper and center panel) we compare the TO position (brightness and colour) in 
the age range 9-14 Gyr for the two sets of isochrones with a c g and a Q . We have transformed the 
isochrones to the observational V-(B-V) plane according to the colours and bolometric corrections 
used by Salaris & Weiss (1998), but the results of this comparison do not depend on the particular 
set of transformations used. As it is evident from the figure, the age differences as derived from 
the TO brightness (or colour) are basically negligible. 

In the lowest panel of Figure 3 we compare A(B — V) as a function of the age, for both 
sets of isochrones. At each metallicity the two curves corresponding to the two calibrations of 
the convection lie parallel for all the relevant age range (the absolute values being different by 
only f«0.02mag). Therefore, the derivative 5(A(B — V))/5t (and the relative ages derived for the 
A{B — V)) is not affected at all when a c g is used instead of a & . 

In conclusion, the main results of this analysis show that the T e g of GC isochrones computed 
employing the MLT formalism and a©, or the a e g calibration as derived from detailed RHD 
models are in good mutual agreement: the maximum deviations in the relevant age range for 
galactic GC amounts at most to a systematic shift by ~ 50±20 K along the RGB. As previously 
discussed, preliminary comparisons (Ludwig et al. 1998) between the adopted grid of 2D RHD 
models and a small sample of 3D ones show only a very small systematic shift of a c g as derived 
from the RHD models by «0.07, which does not affect our results appreciably. 

The next necessary step for finally solving the problem of superadiabatic convection in stellar 
envelopes involves the computations of 3D model grids with up-to-date equation of state and 
frequence-dependent opacity tables to improve especially the photospheric temperature structure. 
Particularly the resolution of the numerical grid in the vertical direction should be improved 
to resolve the extremely sharp temperature jump at the bottom of the photosphere when the 
computations are extended to higher luminosities. A better coverage of the transition region from 
efficient to weak (radiation dominated) convection at high effective temperatures would improve 
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the base to judge between the MLT and the CM formalism. Nevertheless, the already existing 
grid of 2D models indicates that in the T e g -log g- [M/H] region of interest for GC stars there is no 
drastic change in the properties of the envelope convection. Accordingly, the use of the MLT with 
a constant solar-calibrated a© leads only to insignificant errors of at most 0.2 Gyr in the derived 
ages of globular clusters. 

It is important also to remark again that the structure of the superadiabatic convective 
regions is not suitably reproduced either by a & nor by a c g, and that the complete results from 
RHD models have to be employed whenever a detailed description of the properties of this layers 
is needed (e.g., for astro- and helioseismology) . In addition, the RHD models should be analysed 
regarding the effects on the T e g -colour-transformations. 
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Fig. 1. — Upper panel: isochrones in the T e s-log(L / L & ) diagram for [M/H]=— 2.0. Two isochrones 
(t=9, 13 Gyr) computed with a c g are compared to the corresponding isochrones computed with 
a constant a Q . Lower panel: values of a e g along the hydro-calibrated isochrones; the vertical line 
represents the constant a© of the standard isochrones. 

Fig. 2.— As in Figure 1 but for [M/H]=-1.0. 

Fig. 3. — Upper panel: M v values for the TO of isochrones computed by using either a c g or a©, 
[M/H]=-2.0 and -1.0, t=9-14Gyr. Central panel: the same, but for the TO (B-V) colour. Lower 
panel: the same, but for the quantity A(B — V). 
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